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ABSTRACT 

We derive radial profiles of the surface-mass- density for 19 spiral galaxies 
directly from their high-resolution rotation curves. Using the corresponding 
luminosity profiles, we obtain the radial distribution of the mass-to-luminosity 
ratios (M/L) from the inner bulge (~ a few 100 pc) to the outer disk (> 2—10 
kpc) for 11 galaxies (with inclination < 70° in order to reduce the influence 
of the interstellar extinction. The M/Ls in the bulges of two galaxies with 
sufficient resolution, NGC 4527 and NGC 6946, are found to increase steeply 
toward the center at radii ~ 100—500 pc at rates of 15±3 and 7±2 times per 
kpc, respectively. Some other galaxies with fairly high resolution also show signs 
of an increase toward the center. Such an increase may indicate the existence 
of a new component, a "massive dark core", which may be an object linking 
the bulge and a central black hole. Based on radial variations of the M/L, we 
further discuss the variation of the dark-mass fraction in spiral galaxies. 



Subject headings: Galaxies: kinematics and dynamics — Galaxies: photometry 
— Galaxies: spiral — Galaxies: structure 



1. Introduction 



A Radial variation of the mass-to-luminosity ratio in a galaxy (hereafter M/L) is a 
clue to investigate the distributions of visible and dark (invisible) masses (e.g., Kent 1986, 
1987; Persic & Salucci 1988, 1990; Salucci & Frenk 1989; Forbes 1992; Persic et al. 1996; 
Heraudeau & Simien 1997; Rubin et al. 1997). However, most of the current studies have 
been devoted to investigate the dark halo, which is supposed to be closely related to galaxy 
formation and the stability of disks (e.g., Ostriker & Peebles 1973; Athanassoula et al. 
1987). On the other hand, there have been few systematic investigations of dark matter 
within a bulge and an inner disk, for which, however, high-resolution and high-quality 
rotation curves are required as well as corresponding luminosity profiles. 
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Rotation curves (RCs) in disk and outer regions have been derived from position- 
velocity diagrams obtained by optical and Hl-line spectroscopy (Rubin et al. 1980, 1982, 
1985; Bosma 1981a; Mathewson et al. 1992; Persic & Salucci 1995; Mathewson & Ford 
1996; Clemens 1995; Sofue 1996, 1997; Honma & Sofue 1997). Recently, we have obtained 
inner-most RCs by analyzing high-resolution CO, Ha, and [NII]-line spectroscopic data 
(Sofue 1996, 1997; Sofue et al. 1997, 1998, 1999). Fig. 1 shows thus obtained RCs of Sb 
and Sc galaxies (e.g., Sofue 1999). All the Sb galaxies in our sample show very steep rises 
of rotation velocity in the central 100 — 500 pc region. The Sc galaxies also show similar 
nuclear rises, except for a few small-mass galaxies, which show rigid-body growth in the 
RCs. Fig. 1 includes some barred galaxies, which also show similar rotation properties to 
normal galaxies. Recent CCD observations of the Ha-line have also revealed the nuclear 
rises of rotation velocity for many galaxies (Rubin et al. 1997; Sofue et al. 1997, 1998, 
1999; Bertola et al. 1998), in agreement with the CO results. It is also interesting to note 
that the rotation velocities in many galaxies do not appear to decline to zero at the nuclei, 
indicating that the mass density increases steeply toward the center (Bertola et al. 1998; 
Sofue et al. 1999). 

-Fig. 1- 

The fact that almost all the galaxies show steep nuclear rises would imply that the 
high velocities are not due to the non-circular motions in bars. When observing a bar, 
the probability of looking it at a side-on view is fairly higher than that of an end-on view. 
If the influence of the non-circular motion is significant, the RC will indicate lower and 
more gradually rising velocity than the circular velocity, rather close to the pattern speed 
of the bar. Therefore, we may assume that the RCs trace the circular motions in the first 
approximation. If there exist bars, we may still underestimate the circular velocity, which 
would result in an underestimated mass density. 

Radial profiles of the surface-mass density (SMD) and the surface luminosity can be 
used to calculate the M/L directly. As for SMD profiles, Bosma (1981b) used two methods, 
each described by Nordsieck (1973) and Shu et al. (1971), to derive them directly from RCs. 
He compared these methods and found the good agreement between the two SMD profiles. 
However, most of the current investigations have been done on the assumption that the 
visible part of a galaxy consists of a given set of a bulge and a disk, each assumed to have 
a constant M/L. Although this idea is useful to quantify the averaged characteristics in 
each component, it would be necessary to confirm whether the assumption of the constant 
M/L is reliable. Kent (1986) has extensively used the "maximum-disk method" to fit RCs 
calculated on this assumption to observed RCs: He has derived averaged M/Ls in the 
individual components. However, this method is based on the assumption that there is no 
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dark matter in the bulge and disk. Moreover, results of such a fitting process much depend 
on bumps and wiggles in the luminosity profiles (Bosma 1998) as well as on the fitting 
range. 

Forbes (1992) have derived the radial variation of the ratio of the total mass to total 
luminosity involved within a radius, r, which can be called an 'integrated M/U . In this 
paper, we will derive both the SMD and the surface luminosity as differential quantities at a 
given radius, and compare them directly in order to discuss the detailed radial distribution 
of the M/L. 

In section 2, we derive the SMD profiles. We obtain the luminosity profiles in section 
3, describe the obtained M/L profiles in section 4, and discuss the results in section 5. 



2. Radial Profiles of the Surface- Mass Density 

We assume that the 'true' mass distribution in a real disk galaxy will be between two 
extreme cases; spherical and axisymmetric flat-disk distributions. We can calculate the 
SMD, a(R), directly from a RC for these extreme cases. 



2.1. Spherical Mass Distribution 

The mass M(r) inside the radius r is given by 

MM = r -^f, (1) 

where V(r) is the rotation velocity at r. Then the SMD crs(R) at R is calculated by, 

oo 

a s (R) = 2 J p(r)dz, (2) 
o 

1 7 1 dM(r) , 

-dr. (3) 



27r J r\Jr 2 — R 2 dr 



Here, R, r and z are related by r = \^R 2 + z 2 , and the volume mass density p(r) is given by 

P{T) = 47TT 2 ~dr ' (4) 

For a given RC, V(r), eq. (3) can be calculated numerically. Note, nevertheless, an 
underestimation of eq. (3) becomes significant when R ~ -R ma x, the outermost measured 
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radius of the RC, because eq.(3) can be integrated only up to the -R max , instead of oo. 
This edge effect is negligible elsewhere because of the factor (y/r 2 — R?) , by which the 
contribution at r far from R decreased rapidly. In Fig. 2 we show the several results of 
the integration in eq. (3) for the Miyamoto-Nagai potential (Miyamoto & Nagai 1975), by 
taking several different values of -R max . Here, the 'true' value for this potential is given by 
o"d(-R). The figure demonstrates that the edge effect is almost negligible except at radii 
r > 0.9-R ma x, where the SMD is underestimated. 

-Fig. 2- 



2.2. Flat-Disk Mass Distribution 

The SMD for a thin flat-disk, cr D (i?), is derived by solving the Poisson's equation: 



<td(R) 
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(5) 



where K is the complete elliptic integral and becomes very large when x ~ R (Binney & 
Tremaine 1987). 

When we calculate the a-o(R), the following three points must be taken into account. 
(a) First, eq. (5) is subject to the boundary condition, V(0) = V(oo) = 0. The assumption 
V(0) = may not be applied, if there is a black hole (BH) at the galactic center. However, 
a central BH with mass of several to hundreds million M Q dominates the RC only within 
a few pc, and it does not influence the SMD except at the center. On the other hand, 
the spatial resolution in our data is of the order of 100 pc, thus the RCs could hardly 
trace the central BH. So, for our purposes we can safely ignore its existence, (b) Second, 
when calculating the first term on the right hand side of eq. (5), we have only a few data 
points in the central region within a radius comparable to the spatial resolution, where the 
reliability of the calculated V(r) is lower than the outer region, (c) Third, the upper limit 
of the integration of the second term is -R max instead of infinity. Since the RCs are nearly 
flat or declining outward from R = -R max , the second term becomes negative. Hence, we 
overestimate a D (R) at R ~ R max , and, as calculated here, the SMD near the outer edges 
will give an upper limit to the true value. For such galaxies as NGC 1068 and NGC 4569, 
whose rotation velocities appear to increase even outwards of R = R max , the error would be 
still larger. 
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2.3. Verification using the Miyamoto-Nagai Potential 

We, here, examine how the results for the spherical and flat-disk assumptions differ 
from each other as well as from the 'true' SMD by adopting the Miyamoto-Nagai (MN) 
potential (Miyamoto & Nagai 1975), which realistically comprises four components 
denoting a compact nucleus, a bulge, a disk, and a dark halo. Given a set of parameters 
for the MN potential suitable for our Galaxy, we can calculate the 'true' SMD as 
well as the RC. Using this RC, we derive the SMD both for a spherical and flat-disk 
assumptions using the methods as in the previous sections. The MN potential is given 

by $ = Eti GM t [r 2 + {a, + (z 2 + 6 t 2 )°- 5 } 2 ]" ' 5 and ( Ql _ 4 , &!_ 4 , M x _ 4 ) = (0.00, 0.12, 0.05), 
(0.00, 0.75, 0.10), (6.00, 0.50, 1.60), (15.00, 15.00, 3.00) in the units of (kpc, kpc, 10 n M Q ). 

In Fig. 3 we show the calculated true SMD, and the SMDs for a spherical and a 
flat-disk assumptions calculated from the RC. We used the RC only up to R = 20 kpc. Fig. 
3 shows that the spherical case well reproduces the true SMD for the inner region. This is 
not surprising, because a spherical component is dominant within the bulge. On the other 
hand, the flat-disk case mimics the true one in the disk region, which is also reasonable. 
Near the outer edge, the flat-disk case also appears to be a better tracer of the true SMD. 
We stress that the results from the two extreme assumptions, spherical and flat-disk, differ 
at most by a factor of two, and do not differ by more than a factor of 1.5 from the true 
SMD. We may thus safely assume that the true SMD is in between the two extreme cases: 
It is better represented by a spherical case for the inner region, while a flat-disk case is 
better for the disk and outer part. 

-Fig. 3- 

2.4. The Data and Results 

We use the RCs in Sofue (1997) to derive the SMD profiles. We select 18 SA and SAB 
spiral galaxies, for which the spatial resolution (FWHM) is less than 650 pc. Fig. 4 shows 
the calculated <Js(R) and crn(R) with error-bars, respectively. Among the sample, NGC 
4527 has the highest signal-to-noise ratio, which have been obtained from CO-, HI- and 
Ha-line data. The parameters of the sample are listed in Table 1. 

-Tablel- 

-Fig. 4- 

From Fig. 4 the typical difference of SMD values between the two assumptions (sphere 
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and flat-disk) is estimated to be ±0.05 to 0.1 in a logarithmic scale. The error due to the 
error in rotation velocities is about ±0.05. Hence, the effective error of SMD calculation is 
estimated to be ±0.1 to 0.15. 

In Figs. 5 and 6 the SMD profiles of Sb and Sc galaxies are shown together. Although 
some Sc galaxies distinctly have smaller SMD values than the average, the difference in 
SMD between the two galaxy types is fairly small. Namely, Sb and Sc galaxies show similar 
surface-mass distributions. In the disk region at radii 3 to 10 kpc the SMD decreases 
exponentially outward by about 0.3 times/5kpc. In the bulge region (< 2 kpc) of some 
galaxies, the SMD shows a power-law decrease, much steeper than an exponential decrease. 
This may indicate that the surface-mass concentration in the central region is higher than 
the luminosity concentration. The central activities appear to be not directly correlated 
with the mass distribution. Note that the present sample includes galaxies showing Seyfert 
(NGC 1068, NGC 4569), LINER (NGC 3521, NGC 4569), jets (NGC 3079), and black holes 
(the Milky Way). 

-Fig. 5a,b- 

-Fig. 6a,b- 

3. Luminosity Profiles 

3.1. Data Selection 

Luminosity profiles in bulge regions are affected by the interstellar extinction. The 
extinction is higher nearer the center; it is higher for a larger inclination angle. It is also 
higher in shorter wave ranges. In fact, the color-gradient (B — K or V — K) between 
the bulge and disk for galaxies with larger inclinations (i > 60°) is greater than that for 
lower- inclination galaxies (Heraudeau et al. 1996). On the other hand, galaxies with % < 51° 
in de Jong & van der Kruit (1994) do not show significant differences in the color gradient. 

In the present analysis, we have selected the galaxies whose inclinations are less than 
70° if they are observed in NIR-bands and 2-band, and the galaxies with inclinations less 
than 60° if their data are in optical. In the latter case their mean surface luminosities may 
be weakened by 1 mag arcsec~ 2 at most, depending on the inclination. We have, thus, 
selected 11 galaxies out of the 19 RC-selected galaxies. 
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3.2. Description of the Data 

Table 1 lists the parameters for individual galaxies and their references. Since the data 
have been reduced in different ways by different authors, we have checked if their data are 
suitable for the present analysis of Mj L. We briefly describe the data, particularly how the 
surface-luminosity profiles, /i(r), have been obtained. 

(a) Regan & Vogel (1994) have de-projected NGC 598, and calculated yu(r) mag 
arcsec -2 by averaging azimuthally in rings of width ±3". The deviation from axi-symmetry 
is at most about 0.1 mag arcsec -2 including the observational error, which is much smaller 
than the error of the SMD. The effects by spiral arms are small. 

(b) Kodaira et al. (1990) and (d) de Jong & van der Kruit (1994) measured /i(r) along 
the major axes. The errors of fi(r) are as small as ±0.1 — 0.15 mag arcsec -2 . 

(c) Frei et al. (1996) published imaging data of spiral galaxies in FITS format. We 
have used them to derive /i(r) at the same position angles at which the RCs have been 
measured. The errors of /z(r) are found to be smaller than that of SMD. 

(e) Boselli et al. (1997) reduced their data in the same way as Regan & Vogel (1994) to 
derive luminosity profiles. For NGC 4569 they assumed an inclination of 67° and a position 
angle of 23°, which are almost the same as those we adopt here. The errors of /z(r) are less 
than 0.1 mag arcsec -2 between r = 0" and 150". 

(f) Mollenhoff et al. (1995) calculated /x(r) of NGC 4736 along the semi-major axis 
for different position angles and ellipticities. Because of the central weak bar, /i(r) may be 
about 0.5 mag arcsec -2 brighter at r = (10 ± 5) than that at position angle of 108° as 
adopted for the RC. 

We may, thus, consider that the errors in the luminosity profiles caused by the different 
ways is negligible compared to that due to the observations, except for NGC 4736. We 
note, however, that observations in the V-band (NGC 4303, NGC 5055, NGC 6946) may be 
affected by spiral arms, and the errors will amount to about ±0.5 mag arcsec -2 (~ ±0.2 in 
the unit of log(L pc -2 )). 

4. Mass-to-Luminosity Ratios 

To obtain M/L profiles, we first converted the unit of fi(r) from mag arcsec -2 into 
L Q pc -2 , then divided the SMD by {/i(r)(L pc -2 )xcosi} corrected for the inclination. 
We adopt the following solar luminosities in the individual bands: L Q y — 4.83 mag, 
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L &) k = 3.41, L Q K > = 3.45 and L Qji = 4.83. We note, however, that we cannot compare 
the absolute values of M/L profiles directly among the galaxies, because we use different 
filters for different galaxies. Hence, we here focus on relative variations of M/L with radius 
in each individual galaxy rather than compare the absolute M/L values among different 
objects. Fig. 4 shows the thus obtained M/L profiles for the 11 galaxies, and Fig. 7 shows 
the same but for the inner 4 kpc region. As the errors of SMD are much larger than those 
of n(r) for all the galaxies, only the errors in SMD are indicated. 

-Fig. 7- 

4.1. Results for Individual Galaxies 

We discuss the results for the individual galaxies: 

NGC 598: This is a small-mass galaxy which has a mildly rising RC to reach only 
about 100 km s" 1 in the disk region. The M/L increases slowly by 1.7 times from r = 1 to 
3 kpc. In the bulge region it appears to decrease inward rapidly, with large error-bars. It 
also appears to increase inward within 150 pc but it is not conclusive, because the spatial 
resolution is 210 pc. 

NGC 1068: The M/L in the disk region increases slowly at 4.6 times from r = 2 to 6 kpc. 
In the bulge region it decreases more rapidly (0.28 times from r = 2 kpc to 800 pc) than in 
the disk, but within r = 800 pc turns to be nearly flat or somewhat increasing. This is a 
Seyfert type-II galaxy, so the luminosity profile in V-band may be not exact in the central 
region. 

NGC 2903: Both the SMD profile and /i(r) have a remarkable bump at r ~ 2 kpc. Since 
the shape of the RC is also peculiar, we may not neglect the effect of non-axisymmetric 
component such as due to a spiral arm. Interestingly, however, there appears no 
corresponding bump in the M/L. The M/L in the disk region increases slowly by 5.5 times 
from r = 4 to 9 kpc. This Sc galaxy has a compact bulge and we cannot evaluate the M/L 
of the inner bulge for the insufficient spatial resolution. 

NGC 4303: Because of some irregularities caused by spiral arms, the M/L in the disk 
region behaves irregularly, although it increases outward on the average. In the compact 
bulge region the M/L decreases inward steeply from r = 2 kpc to 500 pc. It, then, turns to 
increase within r = 200 pc toward the center. 

NGC 4321: The M/L in the disk increases slowly by 2.0 times from r = 3 to 8 kpc. In 
the bulge region it seems to decrease inward and then increase near the center, though 
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the error-bar is large. Some luminous bumps of about 0.5 mag arcsec -2 are found at r = 
700 pc, 5 kpc and at 11 kpc, which may be caused by spiral arms. The M/L increases 
systematically inward by about 2 times from r = 300 to 150 pc. 

NGC 4527: The M/L in the disk region is almost flat, or slowly increases outward. In 
the bulge region it decreases inward from r = 2 kpc to 500 pc. Then, it turns to increase 
steeply toward the center. 

NGC 4569: The M/L in the disk region increases by 2.4 times from r = 2 to 5 kpc. In the 
bulge region it appears to decrease inward, and then turns to increase toward the center. 
However, the spatial resolution is not sufficient to confirm the central increase. 

NGC 4736: The luminosity profile is obtained only in the bulge region. The M/L in the 
bulge region decreases distinctly inward by 0.38 times from r = 1.8 kpc to 800 pc. The 
Mj L in the inner-most region cannot be evaluated because of the large spatial resolution. 

NGC 5055: The M/L in the disk region increases slowly by 3.0 times from r = 5 to 
15 kpc. In the bulge region its variation is difficult to evaluate because of the insufficient 
resolution as well as due to the large error in the luminosity profile by spiral arms. It looks 
flat in the outer bulge and increases inward in the inner bulge region. 

NGC 5194: The M/L in the disk region is almost flat or rather decreases outward. The 
peculiar behavior of the SMD, and therefore of the M/L, is probably due to a disturbed 
RC by the tidal interaction with the companion galaxy NGC 5195. In the bulge region, the 
M/L is flat or slightly decreases inward to r = 700 pc. Within r = 700 pc, it looks flat or 
slightly increasing, while the spatial resolution is not sufficient to reveal the details. 

NGC 6946: The SMD profile for this galaxy has been calculated with the highest accuracy 
among the 11 galaxies of the present sample (FWHM of 110 pc). The M/L in the disk 
region is flat from r = 1 to 3 kpc and, then, it increases slowly by 2.8 times from r = 3 to 
10 kpc. In the compact bulge region, it shows a remarkable increase toward the center by 
about an order of magnitude from r = 1 kpc to 50 pc. 

4.2. General Characteristics 

From Fig. 4 and Fig. 7 we summarize the general characteristics of the Mj L behavior 
as follows: 

(1) The M/L remains nearly constant in the inner disk. 

(2) It increases gradually outward in the outer disk for all the galaxies, except for NGC 
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5194, which is tidally disturbed by the companion. 

(3) The M/L increases more steeply in the outer halo, indicating the massive dark halo. 

(4) Bulge components show a constant or a slightly smaller value than the disk value. 

(5) For NGC 4527 and NGC 6946 the M/L increases steeply toward the nucleus inside the 
bulge. 

This drastic increase of M/L within r = 500 pc for NGC 4527 and NGC 6946 is of 
particular interest: For NGC 4527 this increase amounts to 4.9 ± 1.0 times/350pc from r = 
500 to 150 pc, and for NGC 6946, 3.5 ± 0.9 times/500pc from r = 600 to 100 pc. Since 
the results are based on optical data, these values may be overestimated because of optical 
extinction. However, it is known that color gradients like V — K and r — K of face-on 
galaxies are of the order of 1 mag arcsec -2 (e.g., de Jong & van der Kruit 1994; Terndrup 
et al. 1994). Therefore, we may conclude that the increase of M/L for NGC 6946, whose 
inclination is only 30°, is significantly greater than that expected from the optical extinction 
effect. For NGC 4527 with % = 60° - 75°, the B - K' color gradient is estimated to be 
about 1.5 mag arcsec -2 (Heraudeau et al. 1996). In the present data, the luminosity profile 
of NGC 4527 is in the Thuan-Gunn i-band, and we may conclude that the central increase 
of M/L for NGC 4527 will be also significant. 

Such a remarkable increase of M/ L toward the center is, however, not detected in the 
other galaxies. For NGC 1068, NGC 2903, NGC 4303, NGC 4569, NGC 4736 and NGC 
5194, the spatial resolution is not sufficient to reveal such an increase even if it did exist. 
For NGC 4321 and NGC 5055 we see some hints of a possible central increase. For the 
small- mass galaxy, NGC 598, we cannot find any increase inward near the center. 

In current studies of the M/L of spiral galaxies, the bulge and disk M/Ls have been 
assumed to be radially constant. However, it does not necessarily mean that the total M/L 
value is constant in the region corresponding to each component: In Fig. 8 we present the 
variations of the M/L for galaxies with inclinations less than 80°. Here, the M/L value is 
normalized to unity at r = 2 kpc. It is now obvious that the total M/L is not constant at 
all within a galaxy, but it varies significantly in the bulge and disk, not only in the massive 
halo. The M/L value often increases outward by 100 times/(10— 30) kpc from the disk to 
the outer halo. Only one exception is NGC 5194, whose peculiar behavior of the M/L 
would be due to the tidal interaction with the companion. 



-Fig. 8- 
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5. Discussion 

5.1. Massive Dark Core 

We have used high-accuracy RC data to derive SMD profiles for 19 galaxies and M/L 
profiles for 11 galaxies. For three galaxies, NGC 598, NGC 4527, and NGC 6946, we 
obtained both high spatial-resolution (r < 200 pc) RCs and luminosity profiles without 
large irregularities due to spiral arms and/or bars. Among these three galaxies, two galaxies, 
NGC 4527 and NGC 6946, are found to show a steep inward increase of M/L in the central 
region, r < 500 pc. Such a steep central increase of M/L may imply that the bulges contain 
an excess of (dark) mass inside a few hundred pc, which we call a "massive dark core" with 
a scale radius of 100 — 200 pc. This yields a mass of M ~ RV 2 /G ~ 10 9 M Q , assuming 
V(R = 100 — 200pc) = 200 km s~\ The massive dark core could be an object linking 
the galactic bulge with a massive black hole in the nuclei (Miyoshi et al. 1995; Genzel et 
al. 1997; Ghez et al. 1998) and/or massive core objects causing a central Keplerian RC 
(Bertola et al. 1998). 

For the Sbc galaxy NGC 4527 we have detected both a dark core and a normal bulge 
within the luminous bulge region. For the Scd galaxy NGC 6946 we found the dark core, 
whereas we could hardly detect a normal bulge component. On the other hand, we could 
see only a normal bulge component in the Scd galaxy NGC 598. Such a difference may 
indicate that the existence of a dark core may not be directly related to the morphology. 
It is also interesting to note that the scale sizes of the dark cores in the two galaxies may 
not be uniquely correlated with the size of luminous bulges. Although dark cores will have 
a crucial implication for the formation, evolution and dynamics of the central bulges of 
galaxies, their origin and universal properties remain open to discussion. 



5.2. Dark-Mass Fraction 

It is interesting to derive the distribution of dark mass fraction (hereafter DMF), -Fdm, 
defined by the local dark-to-total mass density ratio at a given r. If the luminosity profile 
can be corrected for the extinction, we can define the following quantity: 

°"total ^ _ CT* + Q~dark ^ 

where Ctotai, o"*, Cdark, and //* denote the SMD of the total mass, visible mass, dark mass, 
and the surface-luminosity density, respectively. If the M/L of stars on the average at a 
given r, which will be derived from luminosity profiles by multi-band observations in optical 



and infrared bands, is known, then we can define the DMF by 



i'-DM.K.r) — -. — -r-. . {() 

Ctotal 0total//M r J 

This quantity can be used to evaluate the dark-matter distribution. 

Since the absolute calibration of the stellar M/L and the precise correction for the 
extinction are beyond our scope, here we only briefly attempt to examine qualitative trends 
of variation of the DMF. We examine the following two cases: In one case, we assume 
<7*///*(r) to be constant and equal to o"totai// i *( r ) i n the inner disk region; in the other case, 
we take it equal to half of that. We thus derive -Fdm(^) and (Xd ar k(r) for NGC 6946: the 
results are shown in Figs. 9a and 9b, corresponding to the above two cases, respectively. 
Since <7*///*(r) will be smaller in the outer disk region owing to the dominance of younger 
stars of lower metallicity (e.g., de Jong 1996), we may safely consider that the DMF is at 
least 80 — 90 % in the outer disk beyond 8 kpc. In the inner disk region (r ~ 4—5 kpc) the 
DMF appears to be more than 50 %. However, the uncertainty is too large to give a more 
detailed discussion. 

-Fig. 9a,b- 



5.3. Dark Mass v.s. Visible Mass in the Disk Region 

The radial distribution of DMF can be used to derive individual distributions of dark 
mass and visible mass in the disk region. We attempt the following analysis: 

(a) We fit the observed luminosity profiles in a logarithmic scale by a non-linear least 
square fitting by a model, which consists of an exponential bulge and an exponential disk. 
The model is given by 

/Vodei(?") = /^disk(r) + /U bulgc (r), (8) 

where 

/idisk(r) = /x e ,dis k e- L679{r/r — - x) (9) 

and 

/ibui g c(r) = //e.buigee- 1 - 679 ^^-"^- 1 ^ (10) 

We could successfully fit the data by this model for nine galaxies, except for NGC 4736 and 
NGC 5194: NGC 4736 was short of the data in the disk region, and NGC 5194 does not 
have a disk described by the exponential function probably because of the tidal interaction. 
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(b) We define boundary radius, r , at which /U disk //i bulge = 10, and assume that the 
M/L of visible mass, cr*//i*, is constant outside of r . We will focus on the disk region in 
the following. 

(c) We set the M/L value at r = r to that of the total mass, so that the visible 
mass becomes maximum in each galaxy. From the consideration in section 2.3., we adopt 
a flat-disk mass distribution for SMDs. By subtracting the SMDs of visible mass from 
those of total mass, we obtain distributions of dark mass for the nine galaxies. We define, 
"mass-boundary radius", r m b, outside of which the fitted SMDs of visible mass is dominated 
by the dark mass. Note that, because of the assumption of maximum M/L of visible mass, 
we, in fact, derive the maximum mass-boundary radius, r mb max , inside of which dark mass 
may even exceed visible mass. 

In Fig. 10, we show the results of SMD profiles of the total, visible, and dark mass for 
each galaxy. The parameters are listed in Table 2. Two interesting points are found: First, 
dark mass tends to dominate visible mass even in the inner disk region. In Fig. 11, the 
relation between the effective radius of a disk component, r e disk , and r mb max is shown. As 
^mb,max is an upper limit of r m b, r Ci disk is approximately equal to or less than r m b, m ax- This 
trend appears to be independent of morphology, partly because the sample is biased to the 
late-type galaxy. Second, there may be a linear correlation between the two parameters in 
Fig. 11. This trend is caused by the maximum M/L of visible mass in a disk. But the 
M/Ls adopted here range from 0.5 to 4.0 except for NGC 6946, and even if they are set to 
a half, r m b,max will be reduced only by a few kpc. This correlation may indicate that, if any, 
dark mass traces visible mass only up to r = r ej di S k at most, and that the luminous scale of 
galaxies is almost independent of the distribution of a dark halo, because the difference in 
SMD profiles is small (see Fig. 5 and 6). 

-Fig. 10- 

-Fig. 11- 



-Table 2- 
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Figure Captions 

Fig. 1: Top: The rotation curves of 15 Sb (solid lines) and 3 SBc (dashed) galaxies. 
Bottom: Those of 15 Sc (solid lines) and 2 SBc (dashed) galaxies. Note that small-mass 
galaxies, whose maximum disk velocities are 100 to 200 km s _1 , tend to show mild rises. 

Fig. 2: Surface-mass-density profiles calculated from a rotation curve given by 
the Miyamoto-Nagai potential of our Galaxy under an assumption of spherical mass 
distribution. The solid, dashed, and dot-dashed line-profiles represent results for different 
maximum radii, -R max = 20, 18 and 16 kpc, respectively, for the upper limit of the integral 
in eq. (3) in place of oo. 

Fig. 3: Top: The real (solid line) SMD, and spherical (dashed), and flat-disk (dot-dash) 
SMD profiles calculated for the rotation curve of the Miyamoto-Nagai potential. The upper 
limit of the integration is taken as R max = 20 kpc in place of oo. Bottom: The rotation 
curve by the Miyamoto-Nagai model of our Galaxy. 

Fig. 4: Top: (S) Surface-mass- density profile for the spherical assumption. (D) Same 
but for the flat-disk assumption. An bold line denotes the luminosity profile, corrected 
for the inclination. Middle: Mass-to-Luminosity ratio for (S) spherical and (D) flat-disk 
assumptions. Bottom: Rotation curve. The error-bars are caused mainly by those of 
rotation curve. The true values of the surface-mass- density will lie in the shaded area 
between (S) and (D). The horizontal error-bar in the top- left corner denotes the angular 
(spatial) resolution of the surface- mass-density (FWHM), which is larger than that of the 
luminosity profile in all the cases. The interval of the shading lines corresponds to a double 
of the interval of data sampling points. 

Fig. 5a: Top: Surface-mass-density profiles of Sb galaxies calculated for the spherical 
assumption. Because of the edge effect of the integration of eq. (3), the SMD curves decline 
rapidly in the outer 2 to 3 kpc region, where we must not rely on the obtained SMD profiles. 
Bottom: Same but for the inner region at r = to 2 kpc. 

Fig. 5b: Same as Fig. 5a, but for the flat-disk assumption. 

Fig. 6a: Same as Fig. 5a, but for spherical assumption and Sc galaxies. Note that some 
Sc galaxies show smaller SMD values than the total average of our sample. The Scd galaxy 
NGC 598 shows the smallest values in r = ~ 6 kpc. The peculiar behavior of NGC 5194 
at r ~ 8 kpc may be due to the tidal interaction with the companion. 

Fig. 6b: Same as Fig. 6a, but for the flat-disk assumption. 

Fig. 7: M/L profiles at r = to 4 kpc. (a): Calculated for the spherical assumption, (b): 
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For the flat-disk assumption. The meanings of the horizontal and vertical error-bars and 
the shaded area are same as those in Fig. 4. 

Fig. 8: M/L distributions normalized at r = 2 kpc, as obtained by using the SMD and the 
V-band photometric data of Kodaira et al. (1990). 

Fig. 9a: Top: The distribution of the dark mass fraction (DMF) in NGC 6946 for a case 
of <7*///*(r) = 10. The error-bars are due to the difference between the assumptions of 
spherical and flat-disk mass distributions. Note that the DMF is more than 50 % even in 
the disk region. Bottom: The SMD distribution of the dark mass (solid line) and the total 
mass (dashed). 

Fig. 9b: The same as Fig. 9a, but for the <7*///*(r) = 5. Note that more than half of the 
total mass is dominated by dark mass in r > 5 kpc not strongly depending on the cr*///*(r). 

Fig. 10: SMD profiles of the total mass (dot-dash line), visible mass (solid) and dark mass 
(dashed). The narrow lines are based on the two component model, while bold lines are 
directly calculated from the data. 

Fig. 11: The relation between r e>disk and r mb ,max- The straight line shows r mbjmax = r Cidisk . 
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Table 1: (l)Galaxy name, (2) morphological type, (3)inclination (°), (4)distance (Mpc), 
(5),(6)angular (spatial) resolution of the rotation curve, in the units of (") and (kpc), 
(7) photometric band, (8)angular resolution of the luminosity profile ("), and (9)reference 
to the luminosity profile: (a)Regan & Vogel (1994), (b)Kodaira et al. (1990), (c)Frei et al. 
(1996), (d)de Jong & van der Kruit (1994), (e)Boselli et al. (1997), and (f)Mollenhoff et al. 
(1995). Except for NGC 4527, (2)~(6) are taken or calculated from Sofue (1997), while for 
NGC 4527 they are from Sofue et al. (1998). (7) and (8) are taken from their references. 
(5), (6) and (8) are shown as FWHM. 
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Table 2: Fitting parameters. Column 2 and 4 are in the unit of (L pc 2 ) 
, the rest are in the unit of (kpc). 
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